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We study transport dynamics of ultracold cesium atoms in a two-dimensional optical lattice across 
the superfiuid-Mott insulator transition based on in situ imaging. Inducing the phase transition with 
a lattice ramping routine expected to be locally adiabatic, we observe a global mass redistribution 
which requires a very long time to equilibrate, more than 100 times longer than the microscopic time 
scales for on-site interaction and tunneling. When the sample enters the Mott insulator regime, mass 
transport significantly slows down. By employing fast recombination pulses to analyze the occupancy 
distribution, we observe similarly slow-evolving dynamics, and a lower effective temperature at the 
center of the sample. 
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The thorough understanding of atomic interactions in 
optical lattices provides a testing ground to investigate 
hypothetical models widely discussed in condensed mat- 
ter and many-body physics (T], Because of the sim- 
plicity and tunability of the underlying Hamiltonian, re- 
search on optical lattices generates new fronts to perform 
precise, quantitative comparison between theoretical cal- 
culations and measurements. This new class of "precision 
many-body physics" has generated tremendous interest 
in recent years to locate the superfluid (SF) to Mott insu- 
lator (MI) phase boundaries 0-13 , described by the Bose- 
Hubbard model [1,0|) and to characterize Mott and band 
insulators in Fermi gases d, |1] , described by the Fermi- 
Hubbard model [131. Many new, exotic quantum phases 
in optical lattices have also been proposed P] , even in the 
absence of counterparts in condensed matter physics. 

As promising as the precise characterization of quan- 
tum phases is, fundamental assumptions such as the ther- 
mal dynamic equilibrium of the sample should be inves- 
tigated. Since the preparation of quantum gases gener- 
ally involves ramping up the lattice potential, dynamics 
are an inseparable part of all optical lattice experiments. 
Very slow equilibration processes have been reported in 
one-dimensional optical lattices [11] and have been sug- 
gestedby the observation of long-lived repulsively bound 
pairs pJl and doublons [13] in three-dimensional lattices. 
Prospects of non-equilibrium dynamics in optical lat- 
tices have also attracted much interest recently. Mass 
and entropy transport in the optical lattices can provide 
a wealth of information to characterize the underlying 
quantum phases [13, Ell- Dynamic passage across a phase 
transition can lead to the proliferation of topological de- 
fects in the optical lattices [l^ . 

In this letter, we study global dynamics of ultracold 
atomic gases in a monolayer of two-dimensional (2D) op- 
tical lattice. After ramping up the lattice potential, we 
observe both mass transport and statistical distribution 
of atomic occupancy in the lattice. Mass transport is di- 
rectly seen from in situ density profiles, while occupancy 
statistics is probed by inducing loss in sites of three or 
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FIG. 1: (Color online) Averaged absorption images and den- 
sity cross sections of A'^ = 2 x 10'* cesium atoms in a monolayer 
of 2D optical lattice. After ramping the lattice in 150 ms to a 
lattice depth V/ — 13 Er, (a) shows the sample immediately 
after the ramp. In (b), an additional fast recombination pulse 
removes atoms in sites of occupancy three or more, (c) shows 
the average density cross sections of (a) (black circles) and (b) 
(red circles), (d) shows the average density cross sections of 
the samples with additional 800 ms hold time after the ramp, 
without (black circles) and with the recombination loss pulse 
(red circles). Image size is (106 ^m)^ = (200 sites)^ 



more atoms using a fast three-body recombination (see 
Fig. [T|) . Both processes show intriguing behavior at times 
much longer than microscopic time scales for atomic in- 
teraction and tunneling. 

We begin the experiment with a ^'^'^Cs quantum gas 
in a 2D optical trap. Details on the preparation of the 
quantum gas and optical lattice loading procedure can 
be found in Ref . [l^ and Ref . [l^] , respectively. In brief, 
a nearly pure Bose condensate is loaded into a 2D optical 
dipole trap, formed by two orthogonally crossed beams 
on the X — y plane and a one-dimensional vertical op- 
tical lattice of 4 fim spacing which confines the whole 
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sample in a single "pancake" -like lattice site [20j. Us- 
ing microwave tomography, we find ^ 95% of the atoms 
are loaded into a singe pancake trap. The remaining 
^ 5% in the neighboring sites do not contribute to the 
main results reported in this letter. The trap vibra- 
tion frequencies are {ujx,ujy,uJz) —2Trx (11, 13, 1970) Hz, 
and the cloud temperature is T = 11 nK. The ratios 
huJi/kBT = (0.05,0.06,9) indicate the sample is two- 
dimensional. After 2D trap loading, we adjust the atomic 
scattering length a by ramping the magnetic field to a 
designated value, typically, B — 20.7 G where a = 200 as 
and ub is the Bohr radius. At this field, the three-body 
recombination loss rate is at the Efimov minimum [2l|. 

We introduce a 2D optical lattice by slowly turning 
on retro-refiections of the crossed dipole beams which 
add a square lattice potential with lattice spacing d = 
532 nm and a weak contribution to the envelope confine- 
ment characterized by a mean radial frequency ^JuJ^ZJ^ = 
27r(l -f V/%2Er) X 12 Hz, where Er^ ksx 64 nK is the 
recoil energy and V is the lattice depth. Care is taken 
to equalize lattice depths in the x and y directions by 
balancing the lattice vibration frequencies to within 5%. 
Based on the lattice vibration frequency measurement, 
we calculate tunneling t and on-site interaction U nu- 
merically from the band structure in a homogeneous 2D 
lattice. 

We ramp on the lattice depth following V{t) — Vf{l + 

7)/[l + 7e-*(^~^<=)'/^'] [13, preceded by a 30 ms linear 
ramp from to 0.4 Er. The final depth Vf is reached at 
time T = Tc and 7 is chosen such that V{0) = 0.4 Er. 
After the ramp, the sample is held in the lattice for a 
hold time ThoU- The adiabaticity parameter of the ramp 
is given by a — h\t\/t^] slow ramps with a < \ suggests 
that local equilibrium of the system is preserved ITlJlsf . 

We obtain the in situ density profile of the sample 
by absorption imaging normal to the x — y plane. Af- 
ter a hold time Thoid, we first switch the magnetic field 
to B = 17.7 G (a = 40aB) and then turn off the 2D 
lattice 100 /is before the imaging, reducing the on-site 
peak density by a factor of 30, in order to mitigate any 
density dependent loss during the imaging process [22| . 
The atomic density is measured with a spatial resolu- 
tion of 1.3 /im using a long working distance (34 mm) 
commercial microscope objective. The strength and du- 
ration of the imaging pulse are chosen to keep the travel 
distance of the atoms due to the radiation pressure from 
the imaging beam small compared to the depth of focus, 
while maintaining a good signal-to-noise ratio. 

Our first step to study the global dynamics is to watch 
how the density profile comes to equilibrium after a lat- 
tice ramp. Here we employ a ramp which is locally 
adiabatic, but is fast enough to induce detectable mass 
flow. An example is shown in Fig. [5] (a), where after a 
Tc = 20 ms ramp to Vf ^ 10 Er {U/t = 11, a < 0.6) 
the sample of A'' = 2 x 10"* atoms at scattering length 
a = 200 qb gently expands and the peak density slowly 
decreases. This deformation is consistent with the in- 
crease of repulsive atomic interaction in stronger lattice 
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FIG. 2: (Color online) Evolution of the density profile af- 
ter a short lattice ramp. Following a Tc = 20 ms ramp to 
Vf = 10 Er {U/t = 11, a < 0.6), (a) shows the radial den- 
sity profiles measured after hold times of ms (black circles) , 
200 ms (red circles) and 500 ms (blue circles). Inset shows the 
time evolution of A, normalized to the initial value Ao — A(0) 
(solid circles), see Eq. (1), and the single exponential fit. (b) 
shows the profiles measured at Vf = 12 Er (U/t — 20, q < 1). 
The fitted decay times at different depths of Vf are shown in 
(c), where the dashed line marks the critical lattice depth, see 
text. 



confinement. 

To quantify the rate of mass redistribution, we define 
the root-mean-square deviation of a density profile at 
hold time r from equilibrium as 



1/2 



A(r) 



[ni{T) 
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where the sum goes over lattice sites enclosing the sam- 
ple. Hilr) is the mean occupancy of site i at hold time 
T, obtained by averaging over an annular area centered 
on the cloud, containing site i, and with width 1.3 /im 
(20j . fieq^i is the mean occupancy of site i at equilibrium, 
which we obtain from samples that cease to evolve after 
long hold time of Thoid = 500 ^ 800 ms. 

At lattice depths Vf < 10 Er, the sample shows a 
weak breathing mode oscillation in the first 50 ms of hold 
time. After 50 ms, A(r) can be fit by single exponential 
decays with time constants > 100 ms. When the lattice 
depth reaches V/ = 11 Er or higher, the mass flow slows 
down significantly, see Fig. 2 (b) and (c), suggesting that 
the mass transport is suppressed in this regime. The 
crossover behavior near Vf = 11 Er, where U/t w 15, 
is consistent with a recent observation of the suppression 
of superfluidity at U/t — 16 in a 2D optical lattice [sj, 
and quantum Monte Carlo calculations, predicting that 
the SF-MI transition at the tip of the n — 1 Mott lobe 
occurs at U/t « 16.74 in 2D (H]. For Vf > 13 Er, 
even slower dynamics require much longer hold time and 
the slow loss from three-body recombination limits our 
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ability to determine the mass redistribution time scale 

M- 

The slow dynamics throughout the SF-MI regime indi- 
cate that the global thermalization is much slower than 
the microscopic time scales. Indeed, in the range of 
Vf = 6 ^ 13 Eji, tunneling to neighboring sites occurs in 
Tt — h/ zt — 0.6 ^ 3 ms, where z = 4 is the coordination 
number of the 2D square lattice. 

In the second experiment, we investigate the evolution 
of occupancy statistics. For this, we develop a scheme to 
determine the fraction of sites with three or more atoms 
by inducing a fast three-body recombination loss, and 
comparing the density profiles with and without the loss. 
For cesium atoms, extremely fast three-body loss can be 
induced by jumping the magnetic field near an Efimov 
resonance |2l| . where the loss happens much faster than 
atoms tunnel. 

We induce the recombination loss at V/ = 13 by 
jumping the magnetic field to B = 2 G for a duration 
of 1 ms before imaging at S = 17.7 G. The 1/e time of 
the field switching is below 100 /iS. During the switching, 
the magnetic field from the eddy currents is measured by 
microwave spectroscopy and compensated by a controlled 
overshoot of currents in the magnetic coils. At 2 G, the 
three-body loss rate is as high as (20 /is)^^ for 3 atoms 
in one site, much faster than the tunneling rate I/ti — 
(3 ms)~^, and the 1 ms pulse is sufficient to remove all 
the atoms that could participate in the recombination 
process. 

We analyze the dynamics of on-site statistics by first 
ramping the lattice in Tc =300 ms to Vf =13 En {U /t = 
41, a < 0.1) at scattering length a =310 gb and then 
holding the sample for up to 800 ms. Here, the lattice 
ramp is slow enough to ensure negligible subsequent mass 
flow. Density profiles at different hold times, with and 
without the recombination pulse, are shown in Fig. 3 (a- 
c). A larger fractional loss occurs at the central part of 
the sample where the density is higher, as expected; there 
is no apparent loss in the wing. When the sample is held 
longer in the lattice, we observe a smaller fractional loss, 
which suggests that fewer sites are found with three or 
more atoms. 

The evolution of the statistics is best shown in Fig. 3 
(d), where the atom loss. An, is induced by the recom- 
bination pulse after different hold times. A dramatic dif- 
ference is seen near the center with mean occupancy near 
n = 2. Here the loss fraction reaches An/n = 50% imme- 
diately after the ramp, and it slowly declines to merely 
15% after a hold time of Thoid =800 ms. 

To quantitatively model the loss, we assume, starting 
with n atoms in one site, (n modulo 3) atoms remain af- 
ter the pulse. To test this model, we prepare an ideal 2D 
gas by tuning the magnetic field to B — 17.1 G, where 
a « as . We then quickly ramp on the lattice to 30 E^ 
in 10 ms to freeze the on-site occupancy and perform the 
measurement with the loss pulse. For non-interacting 
particles, we expect the occupancy obeys Poisson distri- 
bution. The calculated atom loss, see black solid lines in 



o 
c 

Q. 

o 
o 
o 
c 

(U 



(a) ms 


(b) 200 ms 






•> 

■ 


• 


■ 





600 ms 



50 



Radius r(d) 




0.0 0.5 1.0 1.5 2.0 

Mean occupancy n 



2.5 



FIG. 3: (Color online) Evolution of the on-site statistics in the 
Mott insulator. (iV = 1.6 x lO", = 300 ms, Vf = 13 Er, 
U/t = 41 and a < 0.1). Upper figures show the density 
profiles of the samples held in the final depth Vf = 13 Er for 
Thoid = (a) ms, (b) 200 ms, (c) 600 ms and then imaged with 
(red circles) and without (black circles) the recombination 
pulse. Shaded areas mark the loss fractions, (d) shows the 
loss An versus mean occupancy n measured after different 
hold times (solid circles): ms (black), 200 ms (red), 400 ms 
(green), 600 ms (blue) and 800 ms (magenta). Gray lines are 
the loss derived from an insulator model, see text, assuming 
ksT /U=l (higher curve), 0.5, 0.3, 0.2 and (lower curve). 
The black line, derived from the Poisson distribution, is in 
good agreement with an ideal gas measurement (open circles). 
The inset shows an extended view. 



Fig. 3 (d) and the inset, is in good agreement with our 
measurement. 

Recombination losses measured with interacting sam- 
ples and slow lattice ramps, on the other hand, deviate 
from the Poisson model toward lower values for all mean 
occupancies (see Fig. 3 (d)). This is a general character- 
istic of the strongly interacting gas. 

To gain further insight into the occupancy statistics in 
an insulator, we compare our measurement with an ana- 
lytic model based on a grand canonical ensemble [2^, 113] . 
In deep lattices with t <^U, ksT, the probability for oc- 



cupancy n can be written as P„ = Q e 



> , where 



Hn ~ {U/2)n(n — 1), (3 = l/fc^T, fj. is the local chemical 
potential and Q = g-/3(ff„-Mn) jg the grand partition 
function. The mean occupancy is then n = '^^n ^^'^ 
the loss is modeled as An = n — ^ P„(n mod 3). Cal- 
culations for kBT/U — 1,0.5,0.3,0.2 and are plotted 
in Fig. 3(d). For n < 2.5, all curves show smaller loss 
than does the Poisson model. An insulator at lower tem- 
perature experiences fewer losses because thermal fluctu- 
ation is reduced. At zero temperature, loss only occurs 
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when the occupancy n > 3 is unavoidable, or equiva- 
lently, n > 2. 

Surprisingly, our loss measurements do not follow the 
model with a uniform temperature for up to 800 ms hold 
time. Describing the deviation from a constant temper- 
ature contour by an effective local temperature Te//(r), 
we find the center of the cloud has a lower T^ff ~ 6 nK, 
while for the wing ~20 nK even after 800 ms of hold time. 
This persistent temperature variation across the sample 
suggests that the heat flow is insufficient to establish a 
global thermal equilibrium even after 800 ms hold time. 
This may be aggravated by a large heat capacity of the 
atoms in the wing. 

Both the slow mass and heat flows observed in this 
work raise the issue of describing quantum gases in op- 
tical lattices using a thermodynamic model. We suspect 
that the slow dynamics is partially due to our large sam- 
ple size of (100 sites)^ and the dimensionality of our sys- 
tem, and partially associated with the critical behavior 
of the system. Across the SF-MI transition, the sample 
enters the quantum critical regime, where long equilibra- 
tion times are expected [l6l [28j . Other interesting mech- 
anisms include the long life time of the exited doublon 



|13| . which could slow down statistical redistribution of 
occupancies while supporting mass transport. Moreover, 
the slow recombination loss during the lattice ramp and 
hold time preferentially removes atoms at the center of 
the sample, which could lead to the radial temperature 
gradient observed in our system. 

In summary, we show that the in situ density profiles 
of atoms in a 2D optical lattice provide a viable tool 
for investigating dynamic processes induced by chemi- 
cal potential and temperature imbalance. In both cases, 
we find equilibration times much longer than the micro- 
scopic tunneling time scale. Further investigation into 
these processes and the relevance of our observation to 
the quantum dynamics in the critical regime will be re- 
ported in the future. 
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